De-coupling of Exchange and Persistence Times in Atomistic Models of Glass Formers 
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With molecular dynamics simulations of a fluid mixture of classical particles interacting with 
pair-wise additive Weeks-Chandler- Andersen potentials, we consider the time series of particle dis- 
placements and thereby determine distributions for local persistence times and local exchange times. 
These basic characterizations of glassy dynamics are studied over a range of super-cooled conditions 
and shown to have behaviors, most notably de-coupling, similar to those found in kinetically con- 
strained lattice models of structural glasses. Implications are noted. 



Facilitated dynamics, as encoded in various kinetically 
constrained lattice models (KCMs),[l|, Q implies a gen- 
eral picture of structural glasses in which never-ending 
excitation lines coalesce, branch and percolate through- 
out space-time. 0, 0| Particles cannot move except when 
intersected by these excitation lines. 0, @, 0] As such, 
principal signatures of glassy dynamics can be viewed 
as manifestations of the distribution of times for which a 
particle must wait to first encounter an excitation, and 
the distribution of times between subsequent encounters. 
These are what we have called the distributions of "per- 
sistence" times and "exchange" times, respectively. 0, @, 
7] Definitions and behaviors of excitation lines, persis- 
tence and exchange are most obvious in idealized lattice 
models. Nevertheless, these concepts have proved useful 
for interpreting intermittent and heterogeneous dynamics 



and transport de-coupling 
systems and real systems, 




in continuous force model 
[Tol 11 1 suggesting that 
these concepts are not limited to KCMs. Here, we show 
that indeed, exchange and persistence are well defined in 
terms of classical trajectories of the sort computed from 
molecular dynamics or observed with microscopy, [l2| and 
for a specific continuous force model, the distributions for 
exchange and persistence times behave similarly to those 
we have previously gleaned from KCMs. 

urn 

Figure 1 illustrates our main result. It shows the dis- 
tributions of persistence and exchange times for particle 
displacement events obtained from molecular dynamics 
simulations of an atomistic model. A local dynamical 
event or excitation is defined as a displacement beyond 
a specified cutoff length (see below). As the liquid be- 
comes increasingly super-cooled the two distributions de- 
couple, and the typical persistence time becomes much 
larger than the typical exchange time. This de-coupling 
occurs because excitations or mobility and their asso- 
ciated excitation lines become relatively sparse at low 
temperatures. As such, for low enough temperatures, 
the time extent of space-time regions devoid of excita- 
tion lines (regions that dominate persistence processes) 
is typically very long compared to that for regions bridg- 
ing the width of excitation lines (regions that dominate 
exchange processes). Thus, the origin of this de-coupling 
is the same as that of dynamic heterogeneity. 



This de-coupling is reflected, for example, in observed 
breakdowns [8j in mean-field transport relations like the 
Stokes-Einstein inverse proportionality between diffusion 



constant, D, and structural relaxation time, T a . LJ] In 
particular, when de-coupling occurs, diffusion is much 
faster than would be predicted from D oc l/r a because 
diffusion, being an exchange process, has l/D propor- 
tional to the first moment of the exchange-time distri- 
bution while r a is the first moment of persistence-time 
distribution. [j| 
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FIG. 1: De-coupling of exchange and persistence times in the 
WCA mixture. A local event is defined as a particle moving 
a distance larger than d. We show results for d = 0.5. (Top) 
Distributions of exchange times and of persistence times for 
particle species A and B for various temperatures T. (Bot- 
tom) Ratio of the average persistence time, t p = (t p ), to the 
average exchange time, r x = (t x ), as a function of T. The 
inset shows r p and r x for both species as a function of 1/T. 
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The model [17[ we study is a variation of the binary 
Lennard-Jones mixture of Ref. Il8l, which has been ex- 
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tensively studied as a model super-cooled liquid (see e.g. 
Ref. [Til ) . We modify this system by removing the attrac- 
tive part of the Lennard- Jones interaction by adopting 
the Weeks- Chandler- Andersen (WCA) separation of the 
pair-potential. [20j We consider a mixture of two particle 
species A and B in a cubic simulation box of side length L 
and volume V = L 3 . The potential energy is the sum of 
the pairwise interactions between two particles of species 
fj,,u e {A, B}, 
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if their separation r is less than 2 1 / 6 ct mi/ , and V^„(r) = 
otherwise. Following Ref. [3 we choose oaa = 1, &bb = 
5/6, ctas = (ctaa + ct_bb)/2 and e^A = £bb = sab = 1. 
The particle masses are = 2 and tob = 1. Lengths, 
times and temperatures are reported in units of ctaa, 
{m,Bo\ A j 'e aa) 1 ^ 2 and eaa/^b j respectively. The number 
of particles of species /Lt is iV^ , and the corresponding mole 
fraction is N^/N = N^^Na + N b ). 

Due to the shortness of the interaction range, 2 1 / 6 <r ~ 
1.12cr, each particle interacts with only its nearest neigh- 
bors, significantly reducing the computational overhead 
of calculating forces. We find that this feature makes the 
WCA mixture up to one order of magnitude faster to sim- 
ulate than the corresponding Lennard-Jones mixture. (l7l| 
It helps us study a reasonably large system at signifi- 
cantly super-cooled conditions. 

Figure 2 demonstrates that the WCA mixture has 
the standard phenomenology associated with glass for- 
mation. The figure presents results from two-point 
time-correlation functions for varying temperatures rang- 
ing from the normal liquid regime to the super-cooled 
regime. 2l| Figure 2(top) and Fig. 2(center) show the 
mean-squared displacement (|Ar(i)| 2 ) and the self- 
intermediate scattering function F s (k,t) = (e lk ' Ar (')) for 
k at the peak of the static structure factor k = kg, respec- 
tively, as a function of time. The curves display the char- 
acteristic low temperature features of increasing slowing 
down, plateaus and stretching observed in simulations 
of similar systems. 3, 23 1 Static density correlations are 
stationary over the range of temperatures considered, 
which suggests nothing of thermodynamic significance is 
happening throughout the simulation. The density cor- 
relation functions decay to their uncorrelated values at 
lengths much smaller than the simulation box length, and 
its structure factor shows no anomalies at small wave 
vectors. These facts indicate that the system remains a 
liquid mixture throughout the simulation and does not 
evolve towards a crystallized or phase-separated state. 
Dynamical quantities show typical liquid-like behavior 
when measured over long enough time scales. For in- 
stance, the small wave vector incoherent scattering func- 
tion is fully consistent with diffusive relaxation. This fact 
also supports the conclusion that the system is not evolv- 
ing towards a crystallized or phase-separated state. 17 1 
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FIG. 2: (Top) Mean squared displacement ((Ar) 2 } = 
(|Ar(t)| 2 ) for both particle species A and B as a function 
of time t for various temperatures T. (Center) Incoherent 
scattering function F s (k,t) = (e lk ' Ar O), evaluated at the 
wavevector k — ko of the first peak in the partial struc- 
ture factors, for both species at the same temperatures as 
above. (Bottom) Violation of the Stokes-Einstein relation 
(dashed line) for decreasing temperatures. The inset shows 
the corresponding relaxation times r a , defined by the rela- 
tion F s (ko,r a ) = 0.1, for each species as a function of inverse 
temperature 1/T. 



Figure 2 (bottom) shows the breakdown of the Stokes- 
Einstein relation at low temperatures. At the lowest tem- 
perature we have studied, the de-coupling between the 
diffusion constant and the structural relaxation time is 
about an order of magnitude. This again is similar to 
what is seen in similar systems. [8( A detailed study of 
dynamic heterogeneity in the WCA mixture is given in 
Ref. [I?] and is left to a future publication. 

To obtain the distributions shown in Fig. 1, we monitor 
the time series of events for each particle in the system, 
defining an event as a particle is displaced beyond a cutoff 
length d. Consider particle i. At the initial time t = 0, 
when we start the observation, its position is r,(0). The 
first event time for that particle, t±, is the first time that 
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particle i has moved far enough so that |rj(ti) — r, (0) | = 
d. A second event occurs for that particle at time t = 
ti +<2i when particle i manages to move again a distance 
d, this time from its position at t\; i.e., \vi{tx + t 2 ) — 
r i(ti)\ = d. A third event takes place after a further wait 
t3, and so on. 

This convention establishes a set of waiting times be- 
tween events for the i-th particle, {ti, *2, ts, . . .}. No- 
tice, however, that t\ has an important physical differ- 
ence from i2, tz, The time t\ is the time for the first 

event to take place without condition on when the previ- 
ous event occurred. The times ti, t 3 , . . . are times between 
events. As noted, we call the time t\ a persistence time, 
and the times t 2 , £3, ■ ■ ■ exchange times. [24 1 

Figure 1 (top) shows the distributions of exchange and 
persistence times for the two kinds of particles in the 
WCA mixture, for temperatures ranging from the normal 
liquid regime to the super-cooled one. We choose the 
cutoff length d to be comparable to the particle size, so as 
to probe particle motion relevant to structural relaxation 
and diffusion. As argued in previous work, [3,0] dynamic 
facilitation features emerge only after a suitable coarse- 
graining of short scale motion. The length d plays the 
role of a coarse-graining length. The distributions of Fig. 
1 are for d — 0.5. We have obtained similar results for d 
in the range d = 0.5-1.0. [1| 

At high temperatures, the exchange and persistence 
time distributions coincide. This is the situation when 
the distribution of exchange times is exponential. Q As 
the temperature is decreased the two distributions be- 
come distinct, that for persistence moving towards longer 
times. Fig. 1 (bottom) shows the ratio t p /t x , where 
t p = (t p ) is the average persistence time and r x = (< x ) is 
the average exchange time. At low temperatures the av- 
erage persistence time becomes much larger than the ex- 
change time. This de-coupling mirrors that of the break- 
down of the Stokes-Einstein relation of Fig. 2. 

The de-coupling of exchange and persistence times 
in the WCA mixture is similar to that observed in 
KCMs. @, HH These distributions reflecting the time se- 
ries of individual particle displacement events are non- 
Poissonian. Event times cluster with periods of quies- 
cence followed by periods of high activity. This clustering 
of event times is one way to view dynamic heterogeneity, 
or more precisely, a coexistence in space-time between 
active and inactive dynamical phases. [2 7| Phase separa- 
tion in trajectory space is a distinguishing prediction of 
dynamic facilitation. [27| De-coupling of persistence and 
exchange processes is a consequence of this phase sepa- 
ration. 

By describing dynamics in terms of exchange and per- 
sistence, one is able to picture particle motion as a 
continuous-time random walk. 28| An individual parti- 
cle makes random walk steps at random times drawn 
from the persistence and exchange time distributions. [Bj 
Along with accounting for the breakdown of Stokes- 



Einstein relations, [B| this simple picture also accounts for 
the non-Fickian to Fickian crossover and the shape [29[ 
of the van Hove correlation function in a wide range of 
systems. 11 
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